Light Quark Mass Dependence in Heavy Quarkonium Physics 
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The issue of chiral extrapolations in heavy quarkonium systems is discussed. We show that the light quark 
mass dependence of the properties of heavy quarkonia is not always suppressed. For quarkonia close to an open 
flavor threshold, even a nonanalytic chiral extrapolation is needed. Both these nontrivial facts are demonstrated 
to appear in the decay widths of the hindered Ml transitions between the first radially excited and ground 
state P-wave charmonia. The results at a pion mass of about 500 MeV could deviate from the value at the 
physical pion mass by a factor of two. Our findings show the necessity of performing chiral extrapolations for 
lattice simulations of heavy quarkonium systems. Furthermore, lattice calculations of these transitions would 
also provide a definite answer to the role of coupled-channel effects in heavy quarkonium physics due to virtual 
heavy mesons. 
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Since the discovery of J/ip, the physics of heavy quarko- 
nium is an important tool for testing QCD. Because both the 
charm and bottom quark masses are much larger than the 
nonperturbative scale Aqcd, heavy quarkonia were well de- 
scribed in the framework of potential models. However, in 
recent years this simple picture has been shattered, as quite 
a few charmonium states close to or above the open charm 
thresholds were discovered, and many of their properties are 
not expected from the potential models. For a recent compre- 
hensive review, see Ref. [1]. The spectrum of heavy quarko- 
nium has been intensively studied using lattice simulations, 
using the quenched approximation in the early stages and in 
full QCD in recent years; for example, see Refs. J2H21- While 
most of the calculations focus on the low-lying states, which 
we refer to as the states below open heavy flavor thresholds, 
only a few calculations tackle the problem of higher excited 
states l8UT0l. So far, all calculations of heavy quarkonium in 
full QCD are performed at light quark masses larger than the 
physical values, or equivalently with unphysical pion masses. 
The simulations in Ref. 1 9 ] are only performed at a single pion 
mass M n = 396 MeV. Mixing of charmonia with pairs of 
open charm states are taken into account in Ref. [ 8 ] at three 
different pion masses ranging from 1 GeV down to 280 MeV, 
yet no chiral extrapolation to the physical pion mass was per- 
formed. In addition to the spectrum, there have also been 
lattice simulations of the charmonium (TT] [121 and bottomo- 
nium radiative transitions lfT3l . The quenched approximation 
is used in Ref. ifTTTl . and the calculations of Ref. 1 12] were 
performed at M„ = 485 MeV 



Being bound states of a heavy quark and heavy antiquark, 
heavy quarkonia do not contain any valence light quark. Thus, 
one would naively expect that the light quark mass depen- 
dence of their properties would be suppressed, so that one can 
use a simple linear formula in the light quark masses [remem- 
ber for example that M% ex (m u + rrid) at leading order] 
for chiral extrapolation, as in Refs. (HE) for mass splittings. 
While this is true for low-lying states, a similar simple extrap- 
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FIG. 1 : Schematic diagrams of the creation and annihilation of sea 
quarks in a heavy quarkonium, with solid and dashed lines represent- 
ing heavy quarks and light sea quarks, respectively. 



olation may not be reasonable for higher excited states. The 
purpose of this Letter is to show that dramatic and even nonan- 
alytic dependences in the light quark masses can arise. Hence, 
for the excited states that are close to open flavor thresholds, 
a formula that takes into account the nonanalyticity should be 
utilized for chiral extrapolation. Furthermore, for radiative 
transitions with strong coupled-channel effects, simulations 
at several pion masses are necessary to extract the physical 
results. 

The effects of light quarks in heavy quarkonium systems are 
due to quantum fluctuations of the sea quarks. Sea quark and 
antiquark pairs are created and annihilated in the color singlet 
heavy quarkonium. Low-energy fluctuations can be described 
in the framework of chiral perturbation theory, which is the 
standard tool for chiral extrapolations. The quarkonium states 
can be included as matter fields. Let us focus on the quark 
mass dependence of the quarkonium mass. Two types of sea 
quark fluctuations are schematically depicted in Fig. [TJ type 
(a) is disconnected and suppressed according to the Okuba- 
Zweig-Iizuka rule; type (b) means that the heavy quark (anti- 
quark) and the virtual sea antiquark (quark) can form a color 
singlet state, a heavy meson (antimeson); i.e. a virtual-heavy - 
meson-antimeson pair is created and annihilated after a short 
propagation. There are certainly other contributions, such as 
the doubly Okuba-Zweig-Iizuka suppressed processes that in- 
duce mixing of the heavy quarkonium and a light meson. We 
expect that such contributions are less important, and there- 



fore do not consider them. Type (a) can be parameterized us- 
ing an effective chiral Lagrangian containing unknown low- 
energy constants. The resulting quark mass dependence is an- 
alytic in the light quark masses up to chiral logarithms, see 
Ref. (HI. For instance, denoting the operator annihilating 
a quarkonium field by ip, a possible contribution would be 
proportional to i/j^i/j(x+) in the effective chiral Lagrangian. 
Here, x+ — u^X u ^ + ux u > ( ) is the flavor trace, and 
X = 2Pdiag (m u , m^, m s ) contains the light quark mass ma- 
trix, where B = \(0\qq\0)\/F 2 and F is the pion decay con- 
stant in the chiral limit. The Goldstone boson field <£, which 
contains the pions, eta and kaons in the SU(3) case, are in- 
cluded in u = VU and U = exp (i\/2$/F). At leading 
order 0(M 2 ), this term gives a contribution proportional to 
B(m u + ma + m s ) to the quarkonium mass, and at (D(M%), 
the chiral logarithm M% log M 2 will arise. 

Complexity comes from type (b), which can lead to nonan- 
alyticity as will be shown below. Because the heavy quarko- 
nium states are normally not far from the open flavor thresh- 
olds, the open flavor mesons, at least the ground states, do 
not necessarily decouple in a low-energy effective field the- 
ory (EFT) for heavy quarkonium. In particular, the masses of 
many excited quarkonium states are very close to the thresh- 
olds. In that case, one should consider coupled-channel ef- 
fects due to coupling to the open flavor mesons and an- 
timesons in chiral extrapolation. As an example, let us study 
P-wave charmonium states. They couple to the pseudoscalar 
and vector charmed mesons in an 5-wave with a coupling 
constant g. The self-energy due to coupling to the charmed 
mesons with masses mi and rri2 is expressed in terms of the 
scalar two-point loop function 
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In the rest frame of the charmonium and taking the nonrela- 
tivistic approximation for both propagators, we regularize the 
divergent loop with a three-momentum cutoff A 
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where M denotes the charmonium mass, c = 2/ii2&i2 with 
jiij = vfiiVfij J (rrii + rrij) the reduced mass and 612 = m\ + 
Vfi2 — M. The mass gets renormalized by the real part of the 
self-energy. Writing out the M n -dependence explicitly, we 
get 

M{M„) = Mo (A, M n ) + # 2 raira 2 Re£(M 2 , A, M^), (3) 

where the factor mirri2 is required for correct normalization, 
and 



M (A, M n ) = M (A) + d(X)M 2 + 0(M*) 



(4) 



is the bare mass. (Note that the mass shift due to the virtual 
loops is a scale-dependent quantity, as can be seen in Eq. ([3]), 
and not a physical observable. For phenomenological stud- 
ies of the charmed meson loops in the charmonium spectrum, 



we refer to Refs. A31H61 .) Both the chiral-limit bare mass 

o 

M (A), and the coefficient d(X) depend on the cutoff, since 
the masses of open flavor heavy mesons mi and m2 depend 
on the pion mass. For simplicity, we assume that the heavy 
mesons lie in the same spin multiplet. Up to G(M 2 ), we have 
rrii = mi + hiM 2 /mi [17], where h\ is a dimensionless co- 
efficient of order unity. Notice that 



Vc= W2£i 
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Ml + 0(M*), 
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where fi l2 = wi^/t^i + m^), and 5 = rhi + m^ — M. 
Therefore, for the case with \S\ < R = M 2 /^ 12 , the unitary 
cut in the loop function Eq. ^ cannot be expanded in a poly- 
nomial in M n . Although R c± 20 MeV is small at the physical 
pion mass, it is around 270 MeV for a pion mass of 500 MeV. 
As a result, there will be a cusp due to the nonanalyticity at the 
point M^k) = mi(M 7r ) + 7x12^,,). Nonanalyticity due to 
similar effects in the chiral extrapolation was discussed earlier 
for the A-resonance CUES and the pion form factor l20l . 

We use h\ = 0.44 as determined from the SU(3) 
mass splittings of both the pseudoscalar and vector charmed 
mesons Ell E2- The M n -dependence of M - M (A) for 
the first radially excited P-wave charmonia is plotted in 
Fig. [2] (upper), where A = 0.63 GeV, corresponding to 
M — Mo (A) = for the h' c at the physical pion mass, is 
used; g[ is the coupling of the 2P charmonia to the charmed 
mesons, as defined in [23], and its dimension is mass -1 / 2 . 
Using model values for the masses of Xcc Xci anc ^ ^c at m e 
physical pion mass from Ref. [24], these states are above the 
coupled thresholds at the physical pion mass. Increasing the 
pion mass, the charmed meson masses increase, too. One ex- 
pects the charmonium mass to increase more slowly than the 
charmed meson thresholds. Therefore, for a charmonium with 
a mass slightly higher than the open charm threshold at the 
physical pion mass, the charmonium mass should coincide 
with the threshold at some larger pion mass. After that, the 
open charm mesons cannot go on shell, and a cusp shows up 
because of the end of the unitary cut, as seen in Fig. [2] (upper); 
Xc2 is always below the D*D* threshold, so that there is no 
cusp in the curve for this state. 

For an £-wave charmonium, the nonanalyticity due to cou- 
pling to the pseudoscalar and vector charmed mesons is less 
important, and even invisible. This is because the coupling 
is in a P-wave. Of course, it can couple to a ground state 
charmed meson and an orbitally excited state in an £-wave, 
as considered in Ref. 0. However, the thresholds are far 
from the masses of the IS and 2S charmonia. In this case, 
the square root in Eq. ([5]) can be expanded in a polynomial in 
M 2 . Hence, the M n -dependence of the quarkonium mass is 

given by Eq. (|4|) with redefined Mo (A) and d(X). 

It is instructive to briefly discuss possible hadronic 
molecules with a binding energy much smaller than the pion 
mass. In this case, the bound state can be described by an EFT 
with only contact terms analogous to that for the deuteron; 
for example, see l25l . Then the pion mass dependence of 
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FIG. 2: Upper panel: Pion mass dependence of M — Mo (A) for the 
2P charmonia calculated with A = 0.63 GeV. Lower panel: Pion 
mass dependence of the DD* threshold. 



the mass of the hadronic molecule is dominated by that of 
the masses of the constituents, as argued in Ref. [26]. So, 
in the pion mass range where X(3872) is a DD* bound 
state 1271 , the M n -dependence of its mass should be ap- 
proximately given by that of the threshold, as depicted in 
Fig. [5] (lower) at 0(M 2 ). We will not calculate the devia- 
tion from the threshold due to a small but finite binding en- 
ergy here, but only point out that a loosely bound state can 
easily become unbound by varying the interaction strength. It 
is worth notice that the coupling constant g in Eq. ([3]), which 
controls the strength of the cusp in the M n -dependence of the 
charmonium mass, is also a measure of the hadronic molec- 
ular content of a given state (28J [29) . However, it is obvious 
that a quantitative treatment of the quark mass dependence of 
a hadronic molecule requires a more refined approach than 
that given here. For example, see Refs. l3Qll3TTl . 

The chiral corrections to the heavy quarkonium mass are 
always small compared to the mass in the chiral limit. More 
noticeable is that there exist quantities in heavy quarkonium 
physics whose pion mass dependence is strong. For these 
quantities, chiral extrapolation is mandatory. One can imag- 




FIG. 3: Hadronic loop diagram. Double, solid and wiggly lines de- 
note the charmonia, charmed mesons, and photon, respectively. The 
dashed curves represent the unitary cuts. 



ine that in the mass splittings between two heavy quarkonium 
states, the great bulk of the chiral-limit masses cancels, and 
the chiral corrections are potentially large. However, as seen 
in Eq. Q, it is not possible to give a parameter- free predic- 
tion for its pion mass dependence even at G(M 2 ), since d(X) 
is scale dependent. A prediction can only be made after fit- 
ting the parameters to sufficiently large data. But there are 
indeed quantities whose quark mass dependence is strong and 
can be predicted parameter free. A good example is given by 
the decay widths of the hindered Ml transitions between the 
2P and IP charmonium states. These transitions are shown 
to be dominated by coupled-channel effects l32l based on a 
nonrelativistic effective field theory (NREFT) I23ll33ll34l. 

As shown in Ref. l32l , these transitions are dominated by 
triangle diagrams at the hadronic level, with three intermedi- 
ate charmed or anticharmed mesons (see Fig. [3]). The decay 
amplitude is proportional to the three-point scalar loop func- 
tion. It is convergent, and the nonrelativistic expression reads 
as 12311351 
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where P and q are the momenta of the initial particle and the 
photon, respectively, 623 = rri2 -\- ms -\- q° — M with M the 
mass of the initial particle, J\f = /ii2M23/(167rmi rri2ms), 
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and c is as defined below Eq. ^. For small a, one may expand 
the loop function out 
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It is then clear that two unitary cuts (see Fig. [3]) lead to the 
main contribution of the three-point loop. Corresponding to 
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FIG. 4: Dependence of the widths of various hindered Ml transitions 
on the pion mass. The vertical line denotes the physical pion mass. 



the two cuts, one may define two velocities of the interme- 
diate mesons. The velocity v used in the NREFT power 
counting [23 ] should be understood as the average of these 
two velocities. Following the discussion around Eq. ([5]), we 
can expect cusps in the M n -dependence of the decay widths 
when the mass of the decaying particle coincides with the cou- 
pled threshold. The results for T(x' C 2 ~~ ^ ^c7) and T(h' c — >> 
XcJl) are shown in Fig.|4J where we have neglected the M^- 
dependence of the charmonium masses and used the same 
model value 3908 MeV for M K El as before. 

The cusp in T(h' c — >> XcJl) appears at about M n = 
300 MeV, the same value as in Fig. [2] for h' c . From Fig. [5J 
one finds a strong dependence on the pion mass. The value 
of T(h' c — >> Xco7) at a pion mass of 485 MeV, the same 
value as that used in the lattice simulations for the Ml tran- 
sitions between 5-wave charmonia 1 12], is only about half of 
that at the physical pion mass. This observation highlights 
the necessity of chiral extrapolation of lattice simulations for 
radiative transitions of heavy quarkonia and the necessity of 
small pion masses in the simulations. Otherwise, the un- 
certainty due to unphysical pion mass could be very large. 
Although parameter- free predictions for Ml transitions be- 
tween the 5-wave heavy quarkonia are not possible, as noted 
in Ref. (32lL the charmed meson loops are still expected to 
be crucial (36H37). Thus, the pion mass dependence induced 
by the virtual charmed mesons could introduce a large uncer- 
tainty due to the large pion mass used in lattice simulations. 

One can further study the strange quark mass dependence, 
which translates into the dependence on M K = Bm s , where 
we use the same notation as Ref. [ 38 1. As an example, 
we plot the simultaneous dependence on M^ and Mk of 
r(Xc2 — >■ 7^c) and Y(h' c — )> 7Xco) m Fig. [5] where we used 
M D (* } = M D( *) + 2h 1 M$ c /M D w + 0(M K ). One clearly 

sees the nonanalyticity in both the M n and Mk dependence 
of the latter. 

In conclusion, we have discussed chiral extrapolations in 
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heavy quarkonium physics, especially for the higher excited 
states. These states are close to open flavor thresholds. As 
a result, chiral extrapolation may be nonanalytic. This ob- 
servation is true for any excited hadron whose mass is in the 
neighborhood of an 5-wave-coupled hadronic threshold. For 
such a state, we propose to perform the chiral extrapolation of 
the mass using 



M(M n ) = M + dM* + y/e + fM*, 



(8) 



where M, d, e, and / are parameters to be fit to the lattice 
data. For states far away from any open flavor threshold, e 
will be much larger than fM% so that the square root can be 
expanded, and one may use only the first two terms in the 
above equation up to (D(M%). Furthermore, we find that light 
quark mass dependence is not always suppressed for heavy 
quarkonium systems. As an example, we show that lattice 
results for the decay widths of the hindered Ml transitions 
between P-wave charmonia at a pion mass around 500 MeV 
can deviate by a factor of 2 from the actual values at the phys- 
ical pion mass. Simulations of these transitions would also 
provide a nice test of the NREFT, and would be very use- 
ful in identifying the coupled-channel effects, which might be 
the key to understanding some long-standing puzzles in heavy 
quarkonium systems. If the resulting pion mass dependences 
follow our predictions, they would also allow for extraction 
of the product of coupling constants gig[, which cannot be 
measured directly. 
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